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ABSTRACT 

We examine the Sunyaev-Zeldovich effect in the 7-year WMAP data by cross-correlating it with the Planck 
Early-release Sunyaev-Zeldovich (ESZ) catalogue. Our analysis proceeds in two parts. We first perform a 
stacking analysis in which the filtered WMAP data are averaged at the locations of the 175 Planck clusters. 
We then perform a regression analysis to compare the mean amplitude of the SZ signal, Y^qq, in the WMAP data 
to the corresponding amplitude in the Planck data. The aggregate Planck clusters are detected in the 7-year 
WMAP data with a signal to noise ratio of 16.3. In the regression analysis we find the that the SZ amplitude 
measurements agree to better than 25%: a = L23 ±0.18 for the fit Y^^^ = aFj'oo"*- 

Subject headings: cosmology: the cosmic microwave background — galaxies: clusters: general — cosmology: 
large scale structure of universe — cosmology: observations 



1. INTRODUCTION 

Clusters of galaxies are important objects for studying cos- 
mology and large scale structure formation. In hot clus- 
ters, around 12% of their mass i s in the form of a ho t, ion- 
ized intra-cluster medium (ICM. iMcCarthv et alj|2007 b. The 
ICM can be studied using dire ct X-ray imaging and/or the 
Sunya ev-Zeldovich effect (SZ, ISunvaev & ZeldovichI 119721 
I1980I) . Observations of the latter have grown tremendously in 
recent years thanks to an array of powerful new instruments 
jBirkinshaw & Gull|[T978t lBirkinshaw|[T999t iCarlstrom et all 
I2OO2I) ! 

The thermal SZ effect is a secondary anisotropy in the the 
cosmic microwave background (CMB) radiation, caused by 
CMB photons scattering off free electrons in the hot ICM 
through the inverse Compton effect. The effect boosts the 
photon's energy and thus distorts the CMB spectrum in the 
direction of a cluster, causing a temperature decrement on the 
low frequency side of the CMB peak (specifically, < 217 
GHz), and an increment at high frequencies. 

The SZ effect is especially powerful for studying high- 
redshift galaxy clusters. Since the Compton >'-parameter 
(the integral of the ICM pressure along the line-of-sight) 
does not diminish with distance, and since the CMB has 
a nearly uniform surface brightness, the SZ effect does 
not diminish with increasing redshift. This makes SZ 
surveys especiaUy suitable for fi nding high redshift clusters 
(iPlanck Collaboration VII I''201 1). These cluster samples can 
be used to constrain cosmological models (e.g., from the 
evolution of the mass function) and to probe the physics of 
structure formation (e.g., from cluster scaling relatio ns and 
structural properties. iPlanck Collaboration VIIlll201 ll ). Many 
ongoing experiments are measuring the SZ effect; for in- 
stance, the South Pole Telescope (SPT, Carlstrom et al. 20 11) 
and th e Atacama Cosmology Telescope (ACT, Marriag e et al.i 
1201 Ih are each measuring tens to hundreds of clusters over a 
few hundred square degr ees ([Haiman. Mohr & Holder 2001 



Weller, Battve & Kneissf '2002; Levine, Schulz, & White 
2002; Maiumdar & Mohr 2004; Douspis, Aghanim & Langer 
2006tlShaw. Rudd & Nagai201 1ft . 
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With the mission to precisely measure CMB temperature 
and polarization anisotropy, the Planck satellite was success- 
fully launched by the European Space Agency (ESA) on 14 
May 2009 (Planck Collaborati onll2005h . Planck carries a sci- 
entific payload consisting of 74 detectors sensitive to frequen- 
cies between 25 and 1000 GHz. Planck scans the sky contin- 
uously with an angular resolution between about 30 arcmin 
(FWH M) at the lowest frequency to ab out 4 arcmin at the 
highest (IPlanck Collaboration VIII||201 Ih . Its combination of 
frequency coverage, sensitivity, and angular resolution, en- 
ables it to measure the spatial anisotropy of the CMB with an 
accuracy set by fundamental astrophysical limits. 

Performing an SZ cluster survey over the full sk y was an 
important goal for Planck dAghanim et alj 119971) . and the 
Project produce d its Early Release Compact Source Cata- 
logue (ERCSC, Planck Collaboration VIIll20Tlh in 201 1. The 
ERCSC included a catalogue of 189 SZ clusters detected with 
high-reliability from the first ten months of the Planck survey 
(see IPlanck ERCSC website!) . The Project is ultimately ex- 
pected to release a few thousan d high-reliability SZ clusters 
(IPlanck Collaboration XIIll20Tlh . 

The WMAP mission was launched by NASA on 30 June 
2001 to map the CMB anisotropy over the full sky to multi- 
pole moment / 1000 (angular resolution ^ 0.2°). WMAP 
has precisely measured the cosmological parameters to un- 
precedented accuracy, but its angular resolution and frequency 
coverage (23 to 94 GHz) were not optimized for SZ detection. 
WMAP independently detected a dozen SZ clusters, all of 
which were well-known (iKomatsu et alJl20Tll) . It is therefore 
interesting and important to see if the SZ clusters detected by 
Planck are also detected in the 7-year WMAP data. That is 
the aim of this paper 

The plan of this paper is as follows: in Section|2] we review 
the main features of the SZ effect and discuss the universal 
pressure profile we use to model the ICM pressure. In Section 
[3j we review the Planck ESZ catalogue and the WMAP W- 
band data which forms the basis of our analysis. In Section]?] 
we present the matched filter we apply to the WMAP sky map 
and examine the filtered model profiles for a range of clu ster 
parameters. We present our main results in Section 14.31 and 
some concluding remarks follow. 

Throughout this paper, we adopt a fiducial flat ACDM cos- 
mology with Hubble constant Hq = 70 kms~'Mpc"', and mat- 
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ter density parameter il^ = 0.3. The Hubble parameter at red- 
shift z is H(z) = HqE{z), where E\z) = ^m(i +z)^ + nA. 

2. CLUSTER PRESSURE PROFILE 
2.1. SZ effect 

The thermal SZ effect is a secondary anisotropy in the 
the cosmic microwave background (CMB) radiation, caused 
by CMB photons scattering off free electrons in the hot 
ICM through the inverse Compton effect. The effect boosts 
the photon's energy and thus distorts the CMB spectrum, 
causing a temperature decrement on the low frequency side 
of the CMB peak (specifically, i/ < 217 GHz), and vice 
versa. The induced temperature anisotro py in th e direction 
of a cluster is tSunvaev & Zeldovichtl9'72tlBirkinshawii 19991 
ICarlstrom et al.ll2002h 
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e''+l 
e'' - 1 



(1) 



where = [ri(e^ + l)/(e^ - 1)] -4 captures the frequency de- 
pendence, and 



11 = 



= 1.76 



100 GHz 



(2) 



The dimensionless Comptonization parameter y depends on 
the electron temperature, TJj), and density, «e(r), in the ICM 



y= / «e(r)crT ^al, 



(3) 



where aj is the Thomson cross section, is the Boltzmann 
constant, nif-c^ is the electron rest mass and the integral is 
taken along the line of sight. Since y is positive, the sign of 
determines whether AT is an increment or decrement in the 
CMB temperature. With Tcmb = 2.725 K, we have g,^ >Q for 
v >2n GHz, and vice versa, and gi, — > -2 at low frequencies 

It is convenient to define the integrated SZ signal as F = 
J ydft, where the integration is over the solid angle of the 
cluster. This is equivalent to a volume integral 



1 (7j 



PdV, 



(4) 



where Da is the angular diameter distance to the system and 
P = H ekTe, the electron pressure (PlanckCollabpration VIII 
1201 lb . In the following, we denote the integral performed 
over the sphere of radius R^q^ (SRsoo), as Fgoo (Jsrsoo)- Note 
that Y has units of solid angle, typically arcmin^. 

2.2. Universal profile 

In order to study the SZ effect in the WMAP data, we 
need a model for the pressure profile of the cluster gas, P{r). 
[Arnaud et al] (120101) studied ICM pressure profiles using a 
sample of 33 local (z < 0.2) clusters observed with XMM- 
Newton. The sample spans a mass range of lO'^^Mo < Mgoo < 
IO'^Mq, where Mgoo is the mass enclosed within TJgoo (as- 
sumed to be spherically symmetric). The pressure profiles in 
their sample can be described by a universal profile that is 
scaled with mass- and redshift-dependent factors. The dis- 
persion of the d ata about the scaled profile is less than 30% 
beyond 0.2/^500 (lArnaud et al.ll2010l) . 

^ R500 is defined as the radius within which the density contrast is > 500. 



With X = r/j?5oo, the form of the universal profile given by 
lArnaud et all (120101) is 



P(x)= 1.65 X 10"^£(z)5 



M500 



3 X W^M^h-l 



X p{x) /jyo [ke V cm ■'j , 
where h-jQ = (h/Q.7), ap = 0.12, and 



a' =0.1-(a„ + 0.1) 



l + (x/0.5)3 



(5) 



(6) 



Mere p(x) is tne generaiiz e a JN^' W moaei proposed by 
Naeai. Kravtsov & Vikhliniril (l2007h (see also lAmaud et alJ 
2010) 

P(x) = .... ," ...,«_^^/^ , (7) 



(C500X)T[1+(C50(>X)"]<^-^^/"' 

where Pq, = S.403hjQ is the overall magnitude of the pres- 
sure profile, and C5oo,7,Q:, and /? determine the slope of the 
profile. By fitting this pre ssure template to the simulated 
profile, lAmaud et al.l (120101) found the best parameters are 
C500 = 1 . 1 77 , 7 = . 308 1 , /? = 1 .05 1 , and a = 5 .4905 . 

The only parameter left undetermined in Eq. (|5) is the mass 
parameter M500, which is nearly degenerate with the overall 
pressure normalization term Pq in Eq. (|7]). To illustrate this de- 
generacy, we plot in Fig.[T^ the normalized pressure profiles, 
P(x)/P{Q), for several combinations of M500 and the cluster 
redshift, z- They show very similar normalized profiles. 

The angular profile of the model SZ signal is obtained by 
projecting the 3-dimensional pressure profile P(r) onto the 
plane of the sky and calculating the corresponding tempera- 
ture profile. Following K omatsu et al.. (2011), the projected 
profile, in keV cm"^, is 



Pidid) = ' V (jf + 9^Dl) dl, (8) 

where the Da is the angular diameter distance to redshift z, 
and rout is the truncated radius, which we take to be rout = 
6/?5oo. Beyond this range, the projection is not very sensitive 
to rout because the profile falls off rapidly. 

Given the 2-dimensional pressure profile, the temperature 
profile is 

Arsz(6') CTT 



■.g.^P^m. 
1 nigC 

and the Comptonization parameter y (Eq. ^) is 

y(e)=^P2a{9) 



(9) 



= (8 X 10"^) / dx 
Jo 



xP 



,6'5ooDa/ V^: 



7500 



(10) 



where we have used 9 = r/D^- 

In Fig.[TJ), we plot the normalized Comptonization profile, 
yW/yi^) vs. 9/9500, for the same four cluster parameters as 
in Fig. [T^. Again, the normalized y profiles are very similar 
Since we fit the normalization of the profile to the WMAP 
data, we are primarily concerned with the angular radius of 



WMAP Observations of Planck ESZ Clusters 



3 




each cluster, 6^00, which scales the extent of the profile. As 
discussed below, this parameter may be derived from infor- 
mation provided in the Planck ESZ catalogue. Thus, when 
fitting the WMAP data, we adopt the normalized profile in 
Eq. (doll using the ^500 value predicted by Planck. 

3. DATA DESCRIPTION 

In this section, we introduce the data used in our analysis: 
the Planck Early SZ catalogue, which consists of 189 clus- 
ters; and the WMAP W-band sky map that will be used for 
consistency testing. 

3.1. Planck ESZ catalogue 

Planck is a full-sky CMB survey with nine different fre- 
quency channels: 30, 44, 70, 100, 143, 217, 353, 545, and 
857 GHz. The FWHM angular resolution for these channels 
are 33, 24, 14, 10, 7.1, 5, 5 and 5 arcmin, respectively. 

The ESZ catalogue was obtained from a bUnd, multi- 
frequency search of the Planck-HFl maps, using on an a ll-sky 
extension of the algorithm given in IMehn et al.l (l2006h . The 
team used a matched multi-frequency filter method to enhance 
signal-to-noise in the 189 detected objects; the filter opti- 
mizes detectability by using a linear combination of frequency 
maps to null the CMB signal, and spatial filtering to suppress 
foregrounds and instrument noise ( Planck Collaboration Villi 
1201 II) . In principle, SZ clusters can be also seen in the three 
LFI channels; however, the low-frequency beam size (^25 
arcmin) is generally much larger than a typical cluster size 
of ~5 arcmin, t herefore the LFI channels were n ot included 
in the analysis ( iPlanck Collaboration Villi |201 ll) . The posi- 
tion and radius (5^500) of each cluster profile was varied to 
maximize the signal-to-noise ratio of each detection. The po- 
sition, angular radius, 56*500, and integrated Comptonization 
param eter, F-iffsoo, are tabulated for eac h cluster in the cata- 
logue (IPlanck Collaboration VIIIll20Tllfl . 

Their preliminary analysis yielded 201 high signal-to-noise 
(S/N > 6) candidates, of which 189 were deemed to be of 
high reliability. Of these, 169 were akeady known from X- 
ray and optical surveys, while 20 clusters were newly de- 
tected. The new clusters were subsequentl y confirmed by 
XMM-Newton observations (Planck Collab oration IXl 1201 11) 
and by the AMI survey ( Zwart et al...2OO80 . 

Redshifts are known for 175 of the 189 clusters. Since we 
require redshift information to construct a pressure profile, we 

With the assumption of spherical symmetiy in the model profile, we have 



limit our analysis of WMAP data to this subset of the full cat- 
alogue. The redshift distribution and positions of the Planck 
ESZ clusters are shown in Fig.|2] in the left panel a, the red- 
shift distribution of these clusters is seen to range from 0.01 
to 0.5 with a mean redshift of 0.18. They are distributed uni- 
formly across the sky, outside a \b\ < 14° galaxy cut. 

3.2. WMAP W-band data 

The WMAP satellite produced full-sky maps at five fre- 
quency bands from 23 to 94 GHz, with a FWHM angular res- 
olution from 52.8 to 13.2 arcmin, respectively. These maps 
were used to measure cosmological parameters with unprece- 
dented accuracy. 

Taking into account the combined effects of spectral shape, 
g,y, and diffraction-limited b eam size, the lar gest SZ decre- 
ment occurs around 140 GHz (ICarlstrom et al.t200 2). which is 
higher than the highest WMAP frequency. Within the WMAP 
bands, the 94 GHz W-band map has the most sensitivity to SZ 
signal, owing mostly to its angular resolution (and relatively 
low foreground contamination). In the following analysis, we 
fit SZ profiles to the 7-year W-band sky map. At 94 GHz, the 
SZ decrement is ATsz/T = -I.563'. 

4. FILTERING TECHNIQUE 

To optimally characterize the SZ signal, we need to fil- 
ter the observed maps, which are dominated by primary 
CMB fluctuations, but which also include residual foreground 
signals (including extragalactic point sources) and instru- 
ment noise. In this section we des cribe the choice of filter 
dTegmark & de 01iveira-Costalll998 ), calculate the tempera- 
ture profile of SZ clusters after filtering, and compute the Y- 
parameter We then compare F500 derived from WMAP with 
that derived from Planck. 

4.1. Optimal filter 

SZ clusters are typically unresolved in the WMAP beam, 
so we treat them as the point sources on the sky. In this limit, 
if cluster / has flux 5, at sky position r,, the sky temperature 
Ar(r) may be written as 

AT(r) = cY,SiS(r,ri) + Y,^en,Yem(r), (H) 

/ foil 

where S is the Dirac delta function, and c is the conversion 
factor between flux and temperature, given by 

r2sinh(2)l^ I / he V 10 mK 

V 2^B \kTcMB J MJysr ' 






Fig. 3. — WMAP data used in this analysis. Upper-left - co-added 7-year W-band map. Upp er-r if>ht - WMAP 7-year analysis mask, including point source 
cuts, 82.4% of the sky remains. Lower-left - the filtered W-band temperature map, as per Eq. ( 1151 . Lower-right - the dimensionless y-parameter map, as per 
Eq. (To). 



where 77 is defined in Eq. ^ (see alsolTegmark & Efstathioul 
[1996; Tegmark & de Oliveira-Costa 1998'). Here, the spheri- 
cal harmonic coefficients a^,,, characterize the true CMB tem- 
perature fluctuation. The sky signal convolved with the beam 
response, Bi, is 

+ (13) 

£m 

where fl^°' encodes the temperature due to CMB fluctuations, 
convolved with the beam, and detector noise, not convolved 
with the beam, 

a?^ = agf»B^ + «^„, (14) 



where the beam t ransfer function, Bg, is obtained from 
LAMBDA website. [Note: in this analysis, we use the 
foreground-cleaned maps supplied by the WMAP team, and 
impose the WMAP 7-yr analysis mask (Fig. [3j)); we also 
neglect residual foreground signals in the cleaned, masked 
maps.] Fig. |4] compares the angular power spectrum, Q, di- 
rectly estimated from the map with the spectrum predicted by 
Eq. ( fT4l i using the best-fit ACDM model. The predicted spec- 
trum agrees with the measured spectrum very well. 

To maximize our sensitivity to point sources, we further 
convolve the observed map AT"*"* with an optimal filter Wg, 
so that 



(2£+ 1 \ 
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Fig. 4. — The measured (black) and predicted (red) power spectra, Q, from 
the WMAP 7-year W-band map. The predicted spectrum is based on the best- 
fit ACDM model convolved with the beam, plus detector noise, assuming 
Ni = 0.0187 fj.K^ iHinshaw et al. 2007). 



(15) 



Note that this form implicitly assumes that the beam response 
and optimal filter are both azimuthally symmetric. Treating 
the second line in Eq. (flST l as the noise term, we seek the 
form of Wt that maximizes the cluster signal to noise. In the 
direction of cluster /, the filtered signal is Ar(r,) =ASi, where 



v-^2£-l-l 

: > Bi'We = constant. 



We choose Wi to minimize the ratio 



cr^ = Var 



AT 



(16) 



(17) 



where Cf = B\C^^^+Ni>, and we take C^^^ to be the ACDM 
model power spectrum. 

We minimize Eq. ( fTTj i by adding a Lagrange multiplier to 
the numerator, 

Y.^CTwj^x[Y^^B,wj , (18) 



and minimizing with respect We. We obtain 

Bi Be_ 



B^Lf, 



(19) 



The normalization of We does not affect the signal-to-noise 
ratio of the cluster detection. We plot We in Fig. |5j note that 
We is maximal in the range £ ^ 500-1000. With this filter, 
the smallest variance we can obtain for point source is 



^2^-1-1 Bj 



(20) 



4.2. The filtered cluster profile 

What is the shape of the universal cluster profile after filter- 
ing? Let the unfiltered temperature map due to clusters at 
positions r,-, / = 1 , . . . , A, be 

x(r) = ^Mei), (21) 
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Fig. 5. — The optimal filter (black) for point source detection in the WMAP 
7-year W-band map, Eq. )19t . The CMB signal (red) and detector noise (blue) 
are shown for comparison. 

where /, is the profile of the ith cluster, and 0,- is the angle 
between the /th cluster and r, 

cosG; = ri-r = cos 0; cos 6* -I- sin 6*/ sin (?cos( 0-0,). (22) 

The filtered cluster map may be written as 

x(r) = {W *x){r) 



dn' fi{Q')W(co&Q') 



where 



w(cose') = V WfPKcose'), 

z-^ 47r 



(23) 



(24) 



and COS0' = f ■ r' . In the limit that the SZ clusters can be 



considered point sources, fiiQi) ■ 
reduces to 



: cSiS(ri,r), the filtered map 



x(r): 



v--2£-l-l 
:}2Si22-r—WePe(cosei)- 



4tt 



(25) 



In Fig.|6] we plot selected SZ cluster profiles before and af- 
ter the filtering. For the small cluster case shown in Fig. |6^, 
the filtered profile does not differ appreciably from a filtered 
point source because the cluster is unresolved. For larger clus- 
ter (Fig.|6]3), the filtered profile is noticeably different from a 
filtered point source profile. In either case, we note that 6*500 
lies within the radius where the filtered profile is still posi- 
tive, so this filter should not suppress actual SZ signal. We 
consider this question in more detail in the following section. 

4.3. Results of filtering 

The filtered WMAP 7-year W-band map is shown in the 
lower-left panel of Fig. [3] This map exhibits suppression of 
CMB signal over a range of angular scales. The dimension- 
less y-map obtained from Eq. ( fTOl i is shown in the lower-right 
panel of Fig. [3] We calculate isoo for each cluster and tab- 
ulate the results in column 8 of Tables |2]-2] To evaluate the 
uncertainty in F500, we simulate 1000 sky maps with CMB 
signal and pixel noise. For each map we calculate isoo at each 
cluster position, then compute the standard deviation over the 
ensemble of maps. The resulting error is given in column 9 of 
Tables lam 

Because of beam dilution and detector noise in the WMAP 
data, the signal-to-noise ratio for the detection of individual 
SZ clusters is relatively low; however, the stacked signal is 
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Fig. 6. — The shape of selected cluster profiles before and after filtering. Left - a small cluster with 6500 less than the WMAP W-band beam size, 0.22°; right 
- a marginally resolved cluster. 



TABLE 1 

Aggregate signal-to-noise for stacked SZ clusters 



Method 


WMAP 


Planck 


unweighted, Eq. {26\ 


8.9 


85.2 


weighted, Eq. (271 


16.3 


123.3 



significant. In Table [T] we give the aggregate signal-to-noise 
ratio for two methods of stacking: unweighted and weighted. 
In the first method we evaluate 



(EM) 



1/2 ' 



(26) 



which equally weights all cluster detections. The total signal- 
to-noise ratio for this method is 8.9 for WMAP and 85.2 for 
Planck. In the second method we evaluate 



E 



1/2 



(27) 



which down-weights clusters with low signal-to-noise. This 
gives an aggregate signal-to-noise of 16.3 for WMAP and 
123.3 for Planck. In Fig. [T] we show the aggregate SZ profile 
from the stacked WMAP data and compare this to the pre- 
dicted profile based on stacking the universal pressure profile, 
scaled by the measured Comptonization parameters from the 
Planck catalogue. 

Next, we assess the consistency between the WMAP and 
Planck cluster detections. Figure |8] compares the integrated 
Comptonization parameter, isoo, from each data set. While 
there is quite a bit of scatter, the two measurements are clearly 
correlated. We use linear regression with errors in both axe s 
to fit the slope between the two data sets (see Appendix IaIi. 
Using Eq. iA5\ . we find 

= (1.23±0.18)x yPoo""'' (68% CL), (28) 

which is within 2a of unity. Fig. |8] shows the best-fit regres- 
sion in red. 

While the measured slope is perfectly consistent with unity, 
we briefly consider some potential sources of systematic error 
in our comparison. (1) The amplitude of 5^500 depends on 
how one treats the monopole moment in the WMAP map. In 
our analysis, we subtract the monopole both before and after 
we apply the optimal filter to the map. (2) Finite pixel size can 
introduce noise when identifying pixels within an area that is 
not much larger than the pixel size. We account for this effect 



in our Monte Carlo evaluation of the i'50™'' uncertainty. (3) 
For clusters with 6*500 larger than the WMAP W-band beam 
width of 0.21°, the integral for Fsoo includes some negative 
regions in the filtered profile, as shown in Fig. |6] This will 
suppress the integrated signal in those clusters. The smaller 
value of i'50'™'' measured in the COMA cluster, compared to 
^scM™*^' '-^^^ ^^^^ effect. We note that our value is consis- 
tent with that reported by Komatsu etal] ( 1201 llFI . (4) While 
the Planck team document ed the i r detection of SZ clusters 
jPlanck Collaboration Vil li 120111: iPlanck Collaboration Vll 
I2OI 1; Planck Collaboration IXI I2011h . some processing de- 
tails are not completely clear, such as the precise choice of 
filter. Despite these uncertainties, the best-fit slope between 

and Fsoo™*^ is consistent with unity and suggests that 



.r^wmap 
^500 

such effects are not significant. 

5. CONCLUSION 

The Planck collaboration released its ESZ cluster catalogue 
in Jan 2011. It contained 189 SZ clusters across the full sky, 
175 of which had tabulated redshifts. In this paper, we exam- 
ine these 175 cluster locations in the WMAP data to assess 
consistency with the Planck catalogue. 

We assume the clusters are described by a universal pres- 
sure profile, Eq. (|5}, and project the 3-d profile onto the plane 
of the sky. Given this profile, we can calculate the integrated 
Comptonization parameter, isoo, for any cluster location de- 
termined by Planck. We filter the WMAP 7-year W-band data 
with an optimal filter to suppress primary CMB signal and 
detector noise, then examine the effects of the filter on the 
universal pressure profile. We conclude that the angular ra- 
dius, 6*500, adequately captures the bulk of the SZ signal in the 
filtered data. We estimate the uncertainty in Y^oo using 1000 
Monte Carlo realizations of a ACDM CMB signal and detec- 
tor noise. 

We perform two consistency tests between F50™'' and 

^51)0™''- First, we stack all the cluster data to estimate the total 
signal-to-noise ratio. In the weighted stacking, we obtain an 
aggregate signal-to-noise ratio of 16.3 from the WMAP data, 
which clearly indicates that WMAP detects the most signifi- 
cant clusters seen by Planck. 



^ IKomatsu et all r201 If) do not give an explicit value for 7500, but from the 
measured dimensionless quantity y ~ 7 X 10"^ (two equations after Eq. (71)), 
and the profile width of 10.3' (three lines above Eq. (70) and figure 14), 
we estimate that their measured value of Y^qq for COMA would be ~ 0.01 
arcmin^, which is consistent with our finding within the errors. 
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Fig. 7. — Left - Average SZ profile map obtained from stacl^ing tlie WMAP 7-year filtered W-band data at the locations of the 175 SZ clusters in the Planck 
ESZ catalogue. The color scale ranges from -4 mK (black) to +1 mK (red). Right - The coiTesponding prediction based on stacking the universal pressure profile 
scaled by the PJanci -measured Comptonization parameter, Fjoo- The "Mexican-hat" like ring is the results of applying matched filter. 
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Fig. 8. — Comparison of 175 individual SZ cluster measurements, Ysqo, from WMAP and Planck. The left panel shows the comparison full scale, the right 
panel zooms to a smaller 7500 range. We use the relation 7500 — Jjrsoo/I-S to obtain Y^qq from the Planck catalog. 



Next we compare of and ^500™'' for each cluster In a 
linear regression analysis that accounts for errors in both mea- 
surements, we find the best-fit slope is 1.23 ±0.18, which is 
consistent with unity at the 2a confidence level. We further 
consider some systematic error sources that could lead to a 
slope slightly greater than unity. Our results show that there 
is no fundamental problems with reliability of the ESZ or cal- 
ibration of the SZ amplitu de. Similar conclusion was also 
drawn from 3rd version of IWhitbourn. Shanks & SawangwitI 
(l20Tlh . 

The filtering technique presented here could easily be ex- 
tended to other surveys, such as ACT and SPT We plan to 



revisit this issue when the appropriate data are publicly avail- 
able. 
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APPENDIX 

LIKELIHOOD FUNCTION FOR DATA WITH BOTH X AND Y ERRORS 

In the comparison of yw^^p yPianck^ would like to apply a linear regression model for these two quantities. Since 
both directions have their errors, we need to consider two-dimensional errors when fitting of the slope parameter (see also 
IClutton-Brockll967h . 

Suppose we have one data point (x,y), as plotted in Fig.|9] with measurement errors in both directions (Sx.Sy). The function 
we want to fit is F = f(X), then (X,Y) can be any point on the curve. In practise, (X,Y) should be fairly close to the data point 
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Fig. 9. — Geometric relation between any point on a given function and experimental data point ix.y) with error &y). 

{x^y). We want to measure 

y-f(x) = (y-Y) + (Y-f(x)) 

= (y-Y) + (f(X)-f(x)) 
= (y-Y) + tane(X-x) 
= iy-Y) + f'(x)(X-x). 

The geometric relationship between these quantities have been illustrated in Fig.|9] Therefore, 

Var(y - f(x)) = Var(y - Y) + f'Cx)^ Var((X - x)) 
= Sy^ + f'(xfSx^. 

Therefore, for just one data, we have the likelihood for the function f{x) as 



C{f{x))=- 



[27r {5f- + f'(x)^5x^)\ 



■ exp 



2 I ((5/ + /'(x)2fo2) 



We generalize above equation for a data set with independent data {xii,yk,5xk,5yk), then we have 



a/w)=n 



1 



k=x [2^{5yl+r{xtnxfjY 



- exp 



1 



(yk-f{xk)Y 



2\{5yl + f'{xtnxl 



Therefore, for the case of a linear regression model y = a*x, the likelihood for parameter a becomes 



/:(«)= n 



1 



[2n{5yl + {a5xk)^)\ 



■ exp 



l_ I (yk-a*Xkf 
'2[{5yl + (a6xk)^ 



(Al) 



(A2) 



(A3) 



(A4) 



(A5) 
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TABLE 2 

WMAP OBSERVATIONS OF Planck ESZ CLUSTERS 



Index 


Gal Ion (/) 


Gal lat (b) 


Redshift 


^500 






ywmap 

"soo 


cv^wmap 
•"^500 




(deg) 


(deg) 




(arcmin) 


(arcmin^) 


(arcmin^) 


(arcmin^) 


(arcmin^) 


1 


0.4409 


-41.8351 


0.1651 


6.3071 


0.14£ 


-02 


0.275 


-03 


0.485 


-02 


0.925 


-02 


2 


2.7497 


-56.1827 


0.1411 


7.1114 


0.17£ 


-02 


0.275 


-03 


0.715 


-02 


0.105 


-01 


3 


3.9081 


-59.4171 


0.1510 


7.6784 


0.33£ 


-02 


0.285 


-03 


0.275 


-01 


0.105 


-01 


4 


6.4760 


50.5490 


0.0766 


15.5982 


0.99£ 


-02 


0.845 


-03 


-0.185 


-02 


0.145 


-01 


5 


6.7046 


-35.5407 


0.0894 


10.5796 


0.32£ 


-02 


0.445 


-03 


0.395 


-02 


0.135 


-01 


6 


6.7849 


30.4683 


0.2030 


7.5616 


0.95£; 


-02 


0.405 


-03 


0.255 


-01 


0.105 


-01 


7 


8.3005 


-64.7564 


0.3120 


8.7664 


0.27E 


-02 


0.565 


-03 


0.635 


-02 


0.115 


-01 


8 


8.4486 


-56.3564 


0.1486 


6.7883 


0.13E 


-02 


0.265 


-03 


-0.445 


-02 


0.975 


-02 


9 


8.9362 


-81.2386 


0.3066 


4.6335 


0.23E 


-02 


0.275 


-03 


0.195 


-01 


0.875 


-02 


10 


18.5314 


-25.7228 


0.3171 


4.1424 


0.15£; 


-02 


0.225 


-03 


-0.525 


-02 


0.845 


-02 


11 


21.0918 


33.2560 


0.1514 


9.0695 


0.425 


-02 


0.415 


-03 


0.205 


-01 


0.115 


-01 


12 


29.0054 


44.5625 


0.0353 


21.7792 


0.82E 


-02 


0.125 


-02 


0.565 


-02 


0.155 


-01 


13 


33.4614 


-48.4318 


0.0943 


8.8429 


0.30£ 


-02 


0.345 


-03 


-0.565 


-02 


0.115 


-01 


14 


33.7814 


77.1628 


0.0622 


17.1573 


0.93E 


-02 


0.775 


-03 


0.335 


-01 


0.135 


-01 


15 


36.7219 


14.9232 


0.1525 


7.0561 


0.19£ 


-02 


0.345 


-03 


0.585 


-02 


0.985 


-02 


16 


39.8595 


-39.9889 


0.1760 


5.9780 


0.175 


-02 


0.295 


-03 


0.305 


-02 


0.945 


-02 


17 


42.8257 


56.6172 


0.0723 


12.6087 


0.55£: 


-02 


0.615 


-03 


0.185 


-01 


0.135 


-01 


18 


44.2253 


48.6842 


0.0894 


14.0653 


0.13E 


-01 


0.745 


-03 


0.325 


-01 


0.135 


-01 


19 


46.0815 


27.1816 


0.3890 


3.5573 


0.12E 


-02 


0.195 


-03 


0.875 


-02 


0.765 


-02 


20 


46.5037 


-49.4389 


0.0846 


10.8945 


0.35£ 


-02 


0.425 


-03 


-0.385 


-02 


0.135 


-01 


21 


46.8830 


56.4988 


0.1145 


9.0824 


0.375 


-02 


0.465 


-03 


0.205 


-01 


0.115 


-01 


22 


48.0522 


57.1769 


0.0777 


10.9294 


0.375 


-02 


0.545 


-03 


0.175 


-01 


0.135 


-01 


23 


49.2044 


30.8618 


0.1644 


7.2635 


0.245 


-02 


0.305 


-03 


0.395 


-02 


0.975 


-02 


24 


49.3365 


44.3823 


0.0972 


8.5539 


0.305 


-02 


0.485 


-03 


-0.525 


-03 


0.115 


-01 


25 


49.6670 


-49.5094 


0.0980 


9.5536 


0.215 


-02 


0.375 


-03 


0.185 


-01 


0.125 


-01 


26 


53.4434 


-36.2698 


0.3250 


1.6148 


0.985 


-03 


0.195 


-03 


0.125 


-02 


0.595 


-02 


27 


53.5227 


59.5442 


0.1130 


8.8081 


0.315 


-02 


0.435 


-03 


0.535 


-02 


0.115 


-01 


28 


55.6010 


31.8644 


0.2240 


6.1481 


0.275 


-02 


0.265 


-03 


-0.265 


-02 


0.995 


-02 


29 


55.9755 


-34.8850 


0.1244 


6.7696 


0.205 


-02 


0.285 


-03 


0.655 


-02 


0.985 


-02 


30 


56.8121 


36.3165 


0.0953 


10.6714 


0.325 


-02 


0.405 


-03 


0.755 


-02 


0.125 


-01 


31 


56.9685 


-55.0798 


0.4470 


3.5090 


0.165 


-02 


0.205 


-03 


0.865 


-02 


0.765 


-02 


32 


57.2694 


-45.3577 


0.3970 


4.3122 


0.185 


-02 


0.215 


-03 


0.685 


-02 


0.865 


-02 


33 


57.3362 


88.0113 


0.0231 


40.6731 


0.655 


-01 


0.305 


-02 


0.565 


-02 


0.215 


-01 


34 


57.6138 


34.9421 


0.0802 


10.4783 


0.295 


-02 


0.395 


-03 


0.375 


-02 


0.135 


-01 


35 


57.9289 


27.6442 


0.0757 


11.7297 


0.205 


-02 


0.435 


-03 


0.345 


-01 


0.135 


-01 


36 


58.2818 


18.5938 


0.0650 


12.5308 


0.485- 


-02 


0.505 


-03 


0.585 


-02 


0.135 


-01 


37 


62.4239 


-46.4150 


0.0906 


9.4655 


0.235 


-02 


0.395 


-03 


-0.415 


-02 


0.125 


-01 


38 


62.9261 


43.7096 


0.0299 


27.6348 


0.135 


-01 


0.135 


-02 


0.655 


-02 


0.165 


-01 


39 


67.2317 


67.4641 


0.1712 


7.4168 


0.315 


-02 


0.305 


-03 


-0.305 


-02 


0.105 


-01 


40 


71.6149 


29.7983 


0.1565 


5.4367 


0.135 


-02 


0.205 


-03 


0.625 


-02 


0.925 


-02 


41 


72.6308 


41.4639 


0.2280 


6.2793 


0.475 


-02 


0.275 


-03 


0.125 


-01 


0.975 


-02 


42 


72.8026 


-18.7212 


0.1430 


7.4432 


0.285 


-02 


0.545 


-03 


0.975 


-02 


0.965 


-02 


43 


73.9652 


-27.8216 


0.2329 


6.2334 


0.315 


-02 


0.265 


-03 


0.185 


-01 


0.965 


-02 


44 


77.9099 


-26.6467 


0.1470 


7.2926 


0.225 


-02 


0.305 


-03 


0.125 


-01 


0.105 


-01 


45 


80.3824 


-33.2035 


0.1072 


7.8795 


0.225 


-02 


0.355 


-03 


-0.395 


-02 


0.975 


-02 


46 


80.9954 


-50.9072 


0.2998 


4.5969 


0.155 


-02 


0.255 


-03 


0.135 


-01 


0.895 


-02 


47 


83.2875 


-31.0322 


0.4120 


3.9787 


0.125 


-02 


0.225 


-03 


-0.105 


-01 


0.875 


-02 


48 


85.9999 


26.7107 


0.1790 


11.3240 


0.245 


-02 


0.105 


-02 


0.125 


-01 


0.125 


-01 


49 


86.4555 


15.2999 


0.2600 


4.9172 


0.175 


-02 


0.155 


-03 


0.415 


-03 


0.915 


-02 


50 


92.7308 


73.4614 


0.2279 


5.5586 


0.255 


-02 


0.255 


-03 


0.145 


-01 


0.905 


-02 


51 


93.9197 


34.9077 


0.0809 


11.6216 


0.575 


-02 


0.345 


-03 


0.915 


-02 


0.135 


-01 


52 


94.0187 


27.4256 


0.2990 


8.0500 


0.165 


-02 


0.805 


-03 


0.105 


-01 


0.105 


-01 


53 


96.8523 


52.4668 


0.3179 


4.1455 


0.815 


-03 


0.145 


-03 


0.115 


-02 


0.835 


-02 


54 


97.7396 


38.1199 


0.1709 


6.3915 


0.245 


-02 


0.175 


-03 


0.135 


-01 


0.965 


-02 


55 


98.9502 


24.8610 


0.0928 


8.2040 


0.125 


-02 


0.235 


-03 


0.955 


-03 


0.115 


-01 


56 


106.7311 


-83.2257 


0.2924 


4.2902 


0.205 


-02 


0.245 


-03 


0.115 


-01 


0.835 


-02 


57 


107.1124 


65.3142 


0.2799 


4.9453 


0.175 


-02 


0.205 


-03 


0.215 


-01 


0.885 


-02 


58 


110.9809 


31.7338 


0.0581 


16.6276 


0.135 


-01 


0.505 


-03 


0.145 


-01 


0.145 


-01 


59 


112.4561 


57.0378 


0.0701 


11.2395 


0.295 


-02 


0.355 


-03 


0.155 


-01 


0.135 


-01 


60 


113.8229 


44.3503 


0.2250 


4.6399 


0.655 


-03 


0.135 


-03 


0.655 


-02 


0.875 


-02 
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TABLE 3 

WMAP OBSERVATIONS OF Planck ESZ CLUSTERS 



Index 


Gal Ion (/) 


Gal lat (b) 


Redshift 


^500 








cv^wmap 
•"^500 




(deg) 


(deg) 




(arcmin) 


(arcmin^) 


(arcmin^) 


(arcmin^) 


(arcmin^) 


61 


114.3368 


64.8740 


0.2836 


4.1782 


0.11£ 


-02 


0.18£- 


■03 


0.11£ 


-01 


0.85E 


-02 


62 


115.1624 


-72.0911 


0.0555 


18.8448 


0.12£ 


-01 


0.99£ - 


■03 


0.12£ 


-01 


0.15E 


-01 


63 


1 18.4482 


39.3351 


0.3967 


3.4480 


0.86£ 


-03 


0.14£- 


■03 


0.58£ 


-02 


0.74E 


-02 


64 


118.6016 


28.5586 


0.1780 


5.9658 


0.12£ 


-02 


0.25£ - 


■03 


0.15£ 


-01 


0.96E 


-02 


65 


124.2179 


-36.4859 


0.1971 


6.2310 


0.28£ 


-02 


0.38£- 


■03 


0.66£ 


-02 


0.96E 


-02 


66 


125.5865 


-64.1447 


0.0442 


17.7296 


0.78E 


-02 


0.92£- 


-03 


0.1 8£ 


-01 


0.14E 


-01 


67 


125.7057 


53.8566 


0.3019 


4.1887 


O.lOE 


-02 


0.15£- 


•03 


-0.38E 


-02 


0.87E 


-02 


68 


139.1979 


56.3560 


0.3220 


3.7032 


0.68E 


-03 


0.16£- 


-03 


-0.65E 


-02 


0.79E 


-02 


69 


143.2459 


65.2152 


0.2110 


4.8535 


0.12E 


-02 


0.19B- 


-03 


O.lOE 


-01 


0.90E 


-02 


70 


146.3311 


-15.5913 


0.0172 


39.9005 


0.22E 


-01 


0.33E- 


-02 


-0.28E 


-02 


0.21E 


-01 


71 


149.2421 


54.1894 


0.1369 


7.3566 


0.29E 


-02 


0.26E- 


-03 


O.lOE 


-01 


0.91 E 


-02 


72 


149.7332 


34.6991 


0.1818 


6.9874 


0.33E 


-02 


0.30E- 


-03 


0.78E 


-02 


0.93E 


-02 


73 


159.8592 


-73.4730 


0.2060 


5.6338 


0.29E 


-02 


0.29£- 


-03 


0.87E 


-02 


0.97E 


-02 


74 


161.4447 


26.2337 


0.0381 


17.7570 


0.42E 


-02 


0.68£ - 


-03 


0.22£ 


-01 


0.14E 


-01 


75 


163.7234 


53.5331 


0.1580 


6.4145 


0.17E 


-02 


0.23£- 


-03 


-0.30E 


-02 


0.98E 


-02 


76 


164.1859 


-38.8932 


0.0739 


14.8848 


O.llE 


-01 


O.SIE- 


-03 


0.14E 


-01 


0.13E 


-01 


77 


165.0876 


54.1196 


0.1440 


7.1812 


0.15E 


-02 


0.26E- 


-03 


0.12E 


-01 


0.97E 


-02 


78 


166.1319 


43.3926 


0.2172 


5.6392 


0.21E 


-02 


0.7AE- 


•03 


O.llE 


-01 


0.91E 


-02 


79 


167.6563 


17.6484 


0.1740 


6.3019 


0.25E 


-02 


0.30£- 


•03 


0.89E 


-03 


0.92E 


-02 


80 


171.9474 


^0.6552 


0.2700 


5.9117 


0.34£ 


-02 


0.32£:- 


-03 


0.15E 


-01 


0.94E 


-02 


81 


172.8864 


65.3231 


0.0794 


9.4305 


0.16£ 


-02 


0.33£- 


-03 


0.23E 


-01 


O.llE 


-01 


82 


176.2823 


-35.0506 


0.0347 


25.2705 


0.65£ 


-02 


0.14£- 


-02 


0.14E 


-01 


0.1 6£ 


-01 


83 


180.2409 


21.0459 


0.5460 


3.5482 


0.15£ 


-02 


0.2 1£- 


-03 


0.1 3£ 


-01 


O.llE 


-02 


84 


180.6237 


76.6529 


0.2138 


5.0805 


0.15£ 


-02 


0.21£- 


-03 


-0.9 1£ 


-02 


0.94E 


-02 


85 


182.4440 


-28.2986 


0.0882 


13.0626 


0.92£ 


-02 


0.75£ - 


-03 


0.19£ 


-01 


0.1 3£ 


-01 


86 


182.6361 


55.8245 


0.2060 


5.4790 


O.llE 


-02 


0.235 - 


-03 


-0.38E 


-02 


0.91E 


-02 


87 


186.3949 


37.2555 


0.2820 


5.0003 
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